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Abstract 

Using data from the FOCUS (E831) experiment at Fermilab, we present a new mea- 
surement for the Cabibbo-suppressed decay mode D° — > K + K^tt + -k^ . We measure: 

T(D° -► K+K-it+it-)/T{D q -» K-tt-tt + tt + ) = 0.0295 ±0.0011 ±0.0008. 

An amplitude analysis has been performed in order to determine the resonant 
substructure of this decay mode. The dominant components are the decays 

D° -» Ki(1270)+K-, D° -» Ai(1400) + A~ and D° -» p(770) V(1020). 



1. Introduction 

In recent years there has been a growing interest in multi-body hadronic decays 
of charm mesons. These decays provide a unique tool for investigating the weak 
decay of the charm quark in the environment of low energy strong interactions. 
High statistics data with small backgrounds are still dominated by three-body 
final states, but results on four-body final states are becoming available. 

The picture emerging from many different Dalitz plot analyses reveals a rich 
resonant structure of these decays with a well defined pattern: final states 
that can proceed via simple spectator amplitudes in which the virtual W 
is coupled to a vector (V) or axial-vector (A) meson have large branching 
fractions compared to cases in which the W is coupled to a pseudoscalar (P) 
meson. 
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The decay mode D° — > K + K~tt + tt~ is Cabibbo-suppressed and may be pro- 
duced through two and three-body intermediate resonant states. Following 
the above mentioned pattern, one expects dominant contributions from modes 
having the axial-vector mesons K 1 (D — > AP), namely D° — > K 1 (1270) + K~ 
and D° — > i i fi(1400) + i i f - , both resulting from VU-radiation (spectator) am- 
plitudes. Other tree-level amplitudes lead to intermediate two-body decays 
such as D — > 0(lO2O)p(77O)° (-D — > VV), or even three-body decays like 
L>° -> p(770)°K + K- (D -> VPP). The other possible contribution of the 
D -> W type comes from the mode L>° -> fT* (892)°^* (892)°. In this case, 
however, there is no tree-level amplitude and the two VU-exchange amplitudes 
cancel almost exactly in the SU(3) limit. This is an instance of a mode that 
should proceed primarily through final state interactions. 

Most theoretical predictions of hadronic decay rates are still limited to Cabibbo 
favored modes and two-body decays. There are several predictions for the rate 
of the D° — > p(77O)°0(lO2O) mode [1,2], but we are unaware of any prediction 
for D° -> KtK-. 

In this paper we present a new measurement of the branching ratio T(D° — > 
K + K-it+it-)/T(D -> K-iT-iT + iT + ) using data from the FOCUS experi- 
ment. An amplitude analysis has been performed to determine the D° — > 
K + K^7i + n^ resonant substructure. 

FOCUS, an upgraded version of E687 [3], is a charm photoproduction ex- 
periment which collected data during the 1996-97 fixed target run at Fer- 
milab. Electron and positron beams (typically with 300 GeV endpoint en- 
ergy) obtained from the 800 GeV Tevatron proton beam produce, by means 
of bremsstrahlung, a photon beam which interacts with a segmented BeO tar- 
get [4]. The mean photon energy for reconstructed charm events is ~ 180 GeV. 
A system of three multi-cell threshold Cerenkov counters performs the charged 
particle identification, separating kaons from pions up to a momentum of 
60 GeV/c. Two systems of silicon microvertex detectors are used to track 
particles: the first system consists of 4 planes of microstrips interleaved with 
the experimental target [5] and the second system consists of 12 planes of 
microstrips located downstream of the target. These detectors provide high 
resolution in the transverse plane (approximately 9 /iin), allowing the identi- 
fication and separation of charm primary (production) and secondary (decay) 
vertices. The charged particle momentum is determined by measuring the de- 
flections in two magnets of opposite polarity through five stations of multi-wire 
proportional chambers. 

2. Analysis of the decay mode D° — > K + K ~7r + 7r~ 

The final states are selected using a candidate driven vertex algorithm [3]. A 
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secondary vertex is formed from the four candidate tracks. The momentum 
of the resultant D° candidate is used as a seed track to intersect the other 
reconstructed tracks and to search for a primary vertex. The confidence levels 
of both vertices are required to be greater than 1%. Once the production and 
decay vertices are determined, the distance L between the vertices and its 
error <jl are computed. This is the most important variable for separating 
charm events from non-charm prompt backgrounds. Signal quality is further 
enhanced by cutting on Iso2. This isolation variable requires that all remaining 
tracks not assigned to the primary and secondary vertex have a confidence 
level smaller than the cut to form a vertex with the D candidate daughters. 
To minimize systematic errors on the measurements of the branching ratio, we 
use identical vertex cuts on the signal and normalizing mode, namely L j 
> 9 and Iso2 < 10 %. We also require the primary vertex to be formed with 
at least two reconstructed tracks in addition to the D° seed. 

The only difference in the selection criteria between the D° — > K + K~7t + tt~ 
and D° — > K~tt~7t + tt + decay modes lies in the particle identification cuts. 
The Cerenkov identification cuts used in FOCUS are based on likelihood ratios 
between the various particle identification hypotheses. These likelihoods are 
computed for a given track from the observed firing response (on or off) of all 
the cells that are within the track's (f3 = 1) Cerenkov cone for each of our three 
Cerenkov counters. The product of all firing probabilities for all the cells within 
the three Cerenkov cones produces a x 2 -like variable Wi = —2 ln(Likelihood) 
where % ranges over the electron, pion, kaon, and proton hypotheses [6]. All 
kaon tracks are required to have A K = W w — W K greater than 3 and all 
pion tracks are required to be separated by less than 5 units from the best 
hypothesis, that is picon = W m i n — W w . 

Using the set of selection cuts just described, we obtain the invariant mass dis- 
tribution for K~K + -k~it + shown in Fig. la. Although the Cerenkov cuts con- 
siderably reduce the reflection peak (from D° — > K~7t + 7t~tt + ) to the right of 
the signal peak, there is still a distortion of the background due to this surviv- 
ing contamination. The shape of this reflection peak has been determined by 
generating Monte Carlo D° — > K~7t + 7t~tt + events and reconstructing them as 
K~K + it~it + . The mass plot is fit with a function that includes two Gaussians 
with the same mean but different sigmas to take into account the variation in 
resolution vs. momentum of our spectrometer [3], a second-order polynomial 
for the combinatorial background and a shape for the reflection obtained by 
the Monte Carlo simulation. The amplitude of the reflection peak is a fit pa- 
rameter while its shape is fixed. A log-likelihood fit gives a signal of 2669 ±101 
K~K + ir~7r + events. 

The large statistics K~7t~7t + tt + mass plot is fit with two Gaussians plus a 
second-order polynomial. The fit gives a signal of 131 763 ± 453 K^n^n + n + 
events. 
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The fitted D° masses are in good agreement with the world average [7] and the 
resolutions are in good agreement with those of our Monte Carlo simulation. 



3. Relative Branching Ratio 

The evaluation of relative branching ratios requires yields from the fits to 
be corrected for detector acceptance and efficiency. These differ among the 
various decay modes because of differences in both spectrometer acceptance 
(due to different Q values for the two decay modes) and Cerenkov identification 
efficiency. 

From the Monte Carlo simulations, we compute the relative efficiencies to be: 
n+n+ ) = 0-688 ± 0.006. Using the previous results, we obtain the 
following values for the branching ratio: 

T(D° -»■ K + K-7i + 7i-)/T(D° -»• K-7i-n + 7i + ) = 0.0295 ± 0.0011. 

Our final measurements have been tested by modifying each of the vertex and 
Cerenkov cuts individually. The branching ratio is stable versus several sets 
of cuts as shown in Fig. 2. We varied the confidence level of the secondary 
vertex from 1% to 50%, Iso2 from 10~ 6 to 1, L / a L from 6 to 20, A K from 
1 to 5 and picon from —6 to —2 and consider the tight cut set used to study 
the sub- resonant structure (see below). 

Systematic uncertainties on branching ratio measurements come from differ- 
ent sources. We consider four independent contributions to the systematic 
uncertainty: the split sample component, the fit variant component, the com- 
ponent due to the particular choice of the vertex and Cerenkov cuts (discussed 
previously), and the limited statistics of the Monte Carlo. 

The split sample component takes into account the systematics introduced 
by a residual difference between data and Monte Carlo, due to a possible 
mismatch in the reproduction of the D° momentum and the changing experi- 
mental conditions of the spectrometer during data collection. This component 
has been determined by splitting data into four independent subsamples, ac- 
cording to the D° momentum range (high and low momentum) and the con- 
figuration of the vertex detector, that is, before and after the insertion of an 
upstream silicon system. A technique, employed in FOCUS and modeled after 
the S-f actor method from the Particle Data Group [7] , was used to try to sep- 
arate true systematic variations from statistical fluctuations. The branching 
ratio is evaluated for each of the 4 (= 2 2 ) statistically independent subsamples 
and a scaled error a (that is the errors are boosted when x 2 /(iV — 1) > 1) 
is calculated. The split sample error cr split is defined as the difference between 
the reported statistical error and the scaled error, if the scaled error exceeds 
the statistical error [8]. 
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Source 


Systematic error 


hplit sample 


n no/ 


Fit Variant 




Set of cuts 


1.7% 


MC statistics 


0.9% 


Total systematic error 


2.8% 



Table 1 

Contribution in percent to the systematic uncertainties of the branching ratio 
T{D° -» K-K + 7r-7r + )/r(D° -> if-^-7r+7r+). 

Another possible source of systematic uncertainty is the variant. This com- 
ponent is computed by varying, in a reasonable manner, the fitting conditions 
for the whole data set. In our study we fixed the widths of the Gaussians to the 
values obtained by the Monte Carlo simulation, we changed the background 
parametrization (varying the degree of the polynomial), we removed the re- 
flection peak from the fit function, and we use one Gaussian instead of two. 
Finally the variation of the computed efficiencies, both for D° — > K~ K + -k + it~ 
and the normalizing decay mode, due to the different resonant substructure 
simulated in the Monte Carlo has been taken into account. The BR values 
obtained by these variants are all a priori equally likely, therefore this uncer- 
tainty can be estimated by the r.m.s. of the measurements [8]. 

Analogously to the fit variant, the cut component is estimated using the stan- 
dard deviation of the several sets of cuts shown in Fig. 2. Actually this is an 
overestimate of the cut component because the statistics of the cut samples 
are different. 

Finally, there is a further contribution due to the limited statistics of the Monte 
Carlo simulation used to determine the efficiencies. Adding in quadrature the 
four components, we get the final systematic errors which are summarized in 
Table 1. 

The final result is shown in Table 2 along with a comparison with the previous 
determinations. 



4. Amplitude analysis of D° — > K + K it + it 

A fully coherent amplitude analysis was performed in order to determine the 
resonant substructure of the D° — > K + K~tt + tt~ decay. Previous results on 
this mode were obtained from small samples. E687 [10] did an incoherent 
analysis, while the analysis of E791 [9] quoted only inclusive fractions. 
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Experiment 


r(D°^K+K-TT+w-) 
r(D u ^K-TT-n+n+) 


Events 


r'UUUb (tnis result; 


0.0295 ± U.UU11 ± O.OOOo 


zboy ± 101 


iiwyl [yj 


U.Uolo ± U.UUof ± U.UUoO 


loo ± 10 


E687 [10] 


0.035 ± 0.004 ± 0.002 


244 ± 26 


ARGUS [11] 


0.041 ± 0.007 ± 0.005 


114 ± 20 


CLEO [12] 


0.0314 ± 0.010 ± 0.005 


89 ±29 



Table 2 

Comparison with other experiments. 

Tighter cuts have been applied in the selection of events for the amplitude 
analysis. We require: L / > 10, secondary CL > 0.05, > 4, Iso2 < 10 -5 . 
In addition the secondary vertex must lie outside of the segmented targets to 
reduce contamination due to secondary interactions. Using this set of cuts we 
obtain the invariant mass distribution shown in Fig. lb. The mass plot is fit 
with a function that includes two Gaussians, a second-order polynomial for 
the combinatorial background, and a shape for the reflection peak (similar 
to the fit of Fig. la). A log-likelihood fit finds 1279 ± 48 K~K + ^-^ + signal 
events. 

Figure 3 shows the K~K + , K~ir + , and 7t~tt + projections. We observe clear 
signals of the 0(1020), p(770)° and i^*(892)°. These signals could come from 
decays of the axial- vector mesons -fTi(1270) + and i^i(1400) + , or from decays 
of the type D° -»• VV and D° -»• VPP (V = 0(1020), p(770)°, fsT*(892)°, 
P = ?r ± , i^ ± ). The 7^(1270)+ decay modes that lead to a K~K^^~^ fi- 
nal state are p(770)°K + , K*(892)°tc + , K*(U30)n + and uK + , whereas the 
i^i(1400) + decays only to i^*(892)°7r + . The uK + mode is not considered be- 
cause the branching fraction of uj — > 7r + 7r~ is less than 2%. These components 
are distinguishable by the characteristic angular distributions of the final state 
particles. 

In our model we consider two D — > AP modes, Ki(1270) + K~ and i^i(1400) + fr~ 
We also consider contributions from two D — > VV modes, 0(lO2O)p(77O)° 
and K* (892)° K* (892)°, and from three D -> VPP modes, p(770)°K + K-, 
0(lO2O)vr7T, and K*(B92)°K + n-. Finally, we include the mode / (980)7r+7r- 
(D — > SPP), since the / (980) has a strong coupling to the K + K~ channel. 
The mode fo{980)K + K~ is not considered because there is no phase space 
for it. In principle the mode ao(980)7r + 7r~ could also contribute, but with 
the present statistics one cannot distinguish between the a (980) and / (980). 
Moreover, the dominant component of the a (980) is the r/Ti channel. A good 
description of our data is obtained without the non-resonant channel. 

The formalism used in this amplitude analysis is a straightforward extension 
to four-body decays of the usual Dalitz plot fit technique. The overall signal 
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amplitude is a coherent sum of the ten individual amplitudes, A = J2k c kAk- 
The amplitudes Ak are constructed as a product of form factors, relativistic 
Breit-Wigner functions, and spin amplitudes which account for angular mo- 
mentum conservation. We use the Blatt-Weisskopf damping factors [13], F u as 
form factors (/ is the orbital angular momentum of the decay vertex). For the 
spin amplitudes we use the Lorentz invariant amplitudes [14], which depend 
both on the spin of the resonance(s) and the orbital angular momentum. The 
relativistic Breit-Wigner is 



m 2 — uiq + im r(m) ' 



where 

TV ^ r m ( P *\ 2l+1 F? 
m \PoJ 



In the above equations, m is the two-body invariant mass, m is the resonance 
nominal mass, and p* = p*(m) is the breakup momentum at resonance mass 
m. 

Masses and widths of resonances are taken from the PDG [7], except for the 
p(770)° and the /o(980). The line shape of the p(770)° is taken from Crystal 
Barrel [15], and includes the p — u interference. This significantly improves 
the fit. For the / (980) we used the Flatte formula [16] of a coupled channel 
Breit-Wigner function. The / (980) parameters — g n — 0.20 ± 0.04, g K = 
0.50 ± 0.20 and m = (0.957 ± 0.008)GeV/c 2 — are obtained by a fit to the 
FOCUS D+ -> Dalitz plot, where the / (980)tt + is the dominant 

component. Magnitudes and phases of the _K'i(1270) + decay modes relative to 
the p(770)°K + were determined by the fit, since the limited D° phase space 
will affect the fTi(1270) + decay fractions. 

The fit parameters are the 9 complex coefficients c&. Magnitudes and phases 
are relative to those of the chain D° -> Ki(1270)+FT-, ATi(1270)+ -> p(770)°K + 
The overall signal amplitude is corrected on an event-by-event basis for the ac- 
ceptance, which is nearly constant across the phase space. The finite detector 
resolution causes a smearing of the edges of the five- dimensional phase space. 
This effect is taken into account by multiplying the overall signal distribution 
by a Gaussian factor, g(M) (M being the K + 'K~tt + tt~ mass). The normalized 
signal probability distribution is, thus, 

P S {4>) = ^e{<t>) P {<t>)9{M) |£c fc A fc (0) 2 
with being the coordinates of an event in the five dimensional phase space, 
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the acceptance function, and p(<f>) the phase space density. 

We consider four types of background events: random combinations of a 0(1020) 
and a it~it + pair, a if* (892)° plus a K~ir + pair, a p(770)° plus a K + K~ pair, 
and random combinations of K + K ~7r + 7r~. The relative fractions of these 
backgrounds were determined from a fit to the data on the side bands of 
the K + K~7t + tt~ mass spectrum. This fit yields 69% for random combina- 
tions of K+K-iT+iT-. The fractions of the 0(lO2O)7r+7r-, K* (892)° K~ >k + , and 
p(770)°K + K~ backgrounds are 17%, 5% and 9%, respectively. 

We assume the random K + K~tt + tt~ combinations to be uniformly distributed 
in phase space, whereas for the other backgrounds we assume Breit-Wigners 
with no form factors and no angular distribution. The overall background 
distribution is a weighted, incoherent sum of the four components described 
above. The relative background fractions, bk, are fixed in the fit. The overall 
background distribution is also corrected for the acceptance (assumed to be the 
same as for the signal events) on an event-by-event basis and multiplied by an 
exponential function b(M), accounting for the K + K^n + 7i^mass distribution 
of the background. The normalized background probability distribution is 

P B {<t>) = ^e(0p(06(M)5>fl fc (0. 

An unbinned maximum likelihood fit was performed, minimizing the quantity 
w = — 21n(£). The likelihood function, £, is 

L= II [Ps(P) + Pb(P)\ ■ 

events 

Decay fractions are obtained from the coefficients Ck, determined by the fit, 
and after integrating the overall signal amplitude over the phase space [17]. 
Errors on the fractions include errors on both magnitudes and phases, and are 
computed using the full covariance matrix. 

The result from the best fit is shown in Table 3. The dominant component 
is indeed the mode D° -»• K^(1270 + )K~, with a fraction of (33 ± 6 ± 4)%. 
Integrating the amplitude squared for the summed D — >• AP components 
gives nearly 55% of the total decay rate. This result is consistent with other 
D decays, for instance, the case of D° — > K~ir + ir~ir + , which is dominated by 
the mode ai(1260) + K~ with a fraction of about 50%. 

The second dominant component is the D — > VV mode: a large contribution 
of the 0(lO2O)p(77O)° channel, with a fraction of 29%, and a small, but signif- 
icant, contribution from K* (892)° K* (892)° . The 0(lO2O)p(77O)° channel can 
proceed through an internal W radiation amplitude, while the K* (892)° K* (892) 
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results from final state interactions. Using the value of T(D° — >■ K + K 7t + tt ) /T(D° 
K~it-it + it + ) we measured, we obtain B(D° -> 0(lO2O)p(77O)°) = (1.2±0.1) x 
10~ 3 . This is a factor of five higher than the prediction of Bedaque, Das and 
Mathur [2], (B(D° -> 0(lO2O)p(77O)°) = 2.2 x 10~ 4 ). The model of Bauer, 
Steck and Wirbel [1] predicts B(D° -> 0(lO2O)p(77O)°) = 4.5 x 10~ 3 , which 
is a factor of 3.5 higher than our value. 

Finally, the remaining fraction comes from the D — > VPP and D — > S*PP 
decays. Altogether, these modes account for nearly 30% of the total decay 
rate. A similar result was found in our analysis of the D° — > K + K~K~ir + 
decay [17], showing the importance of the three-body channels in four-body 
decays of D mesons. In addition to the three D — > VPP modes listed above, we 
have also included the mode K*(1400) + K~ , but its contribution is negligible. 

The K + K~ spectrum has an interesting feature. Near the threshold it is dom- 
inated by the 0(1020) peak, but the 0(1020) line shape is distorted by the 
presence of the /o(980). The 0(1020) and / (980) cannot be distinguished 
with purely a mass cut on K + K~ invariant mass. The fairly large fraction 
(15%) of the D° -> /o(980)7r+7r- mode shows that the / (980) contribution 
cannot be neglected. The distinction between 0(1020) and / (980) components 
requires a full angular analysis. 

In four body decays the phase space is 5-dimensional, so we can only look 
at projections. In Figure 4, the K + K~, 7r + 7r~, and K~n + invariant mass 
projections of events used in the amplitude analysis are superimposed on the 
fit result, with the background projections shown in the shaded histograms. 
The fit result can also be displayed in two dimensional projections, shown in 
Figure 5, and in distributions of the cosine of helicity and acoplanarity angles, 
shown in Figure 6. In the case of the D° — > 0(lO2O)p(77O)° mode, the helicity 
angle is defined as the angle between the K~ and the direction of the 
recoiling p(770)° (0(1020)) in the 0(1020) (p(770)°) rest frame. In the case 
of D° — > K* (892)°K* (892)°, helicity angles are defined as the angle between 
the K and the recoiling i^*(892)° in each i^*(892)° frame. The acoplanarity 
angle is the angle between the 0(1020) and p(770)° decay planes, measured in 
the D rest frame. The bump in the central part of the p(770)° helicity angle 
distribution is due to the decay chain D° -> ir 1 (1270)+K'-, fTi(1270)+ -> 
p(770)°^+. 

The invariants used to define the kinematics of the D° — * K + 'K~7t + tt~ decay 
are the K~K + , ir + ir~, K~ir + , K + ix~ and K~n~ masses squared. Due to the 
limited statistics we have integrated over the latter invariant and divided the 
remaining four into three bins each, yielding a total of 81 cells. From these 
a total of 45 cells lie within the phase space of this decay. Data and Monte 
Carlo samples were divided into these 45 cells. The fit has 18 free parameters, 
so the number of degrees of freedom is 45 — 18 = 27. We obtain a x 2 °f 40.4, 
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A/T n o"n l \ 1 1 n p 

IV-Ldglll L U.U.C 


T-* h pi qp 


± i uiuii i /u j 


Ai(lz7UJ A , Ai ^ /9( / /Uj A T 


1 (nxeci) 


U (nxea j 


1 O I /? I o 

18 ± ± 3 


Ai(lz7U) A , Ai — > Aq(14o(JJ7T^ 


r» T7 i c\ no i n n/^ 

0.27 ± 0.U8 ± U.Ub 


354 ± 19 ± 19 


o i 1 i n 

2 ± 1 ± 


Ki{L27[)y K , Ai — > A (892) 7r^ 


U.y4 ± 0.1b ± 0.13 


1 O | 1 Q | 1 FT 

12 ± 12 ± 15 


lb ± 4 ± 5 


Ki{L27{jy A , (all modes) 






OO I /? I /I 

33 ± b ± 4 
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25y ± 11 ± 13 
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oo i o i in 

28 ± 13 ± 1U 
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3 ± 2 ± 1 
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Art i 1 i i o 

4y ± ii ± iz 
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zy ± z ± i 
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U.oo ± U.lz ± U.lo 


z/8 ± ZD ± zU 


O _1_ O _1_ O 

z ± z ± z 


(/>(1020)7r+vr- 


0.3U ± U.U6 ± U.06 


163 ± 16 ± 15 


1 ± 1 ± 


A*(892)°A+7r- 


o.83 ± o.uy ± 0.1U 


234 ± 10 ± 11 


11 ± 2 ± 1 


/ (980)7r+vr- 


o.yi ± 0.13 ± 0.05 


240 ± 11 ± 17 


15 ± 3 ± 2 



Table 3 



Results from the best fit. The second error on the fractions, magnitudes and phases 
is systematic. The fraction for the mode Ai(1270) + A~ (fourth row) includes all 
three decay modes of the A'i(1270) + , added coherently. 

and from this value the estimated confidence level of the fit was 4.7%. 

We consider three different sources of systematic uncertainties for the ampli- 
tudes analysis: split sample, using the same subsamples described previously, 
fit variant, varying the fitting conditions of the whole data set and the compo- 
nent due to the particular choice of the vertex and Cerenkov cuts. The most 
important contributions from fit variant systematic errors are: parameteriza- 
tion of the / (980) line shape, uncertainty in the relative background fractions, 
and uncertainty in the masses and widths of the resonances. We have also con- 
sidered the effect of ignoring form factors and using a flat acceptance. The cut 
component systematic errors were estimated using the standard deviation of 
several different sets of vertex/particle identification cuts, as we did in the 
branching ratio measurement. Systematic errors on phases and fractions were 
obtained adding in quadrature the split sample, fit variant, and cut component 
errors. 



5. Conclusions 

Using data from the FOCUS (E831) experiment at Fermilab, we studied the 
Cabibbo-suppressed decay mode D° — > K + K~it + it~ . 
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A comparison with the two previous determinations of the relative branch- 
ing ratio T(D° — > K + K~7r + 7r~) /F(D — >■ K~7t~7t + tt + ) shows an impressive 
improvement in the accuracy of this measurement. 

A coherent amplitude analysis of K + K~ty + 'k~ final states was performed, 
showing that the dominant contribution comes from D° — > AP modes, cor- 
responding to nearly 55% of the D° — > K + K~tt + tt~ decay rate. We also 
measured an important contribution from D° decaying to two vector mesons, 
corresponding to nearly 30% of the decay rate. The remaining fraction of the 
D° — > K + K~n + ix~ decay rate comes from decays of the type D° — > VPP 
and D° —> SPP, with important contributions from K*(892)°K~tt + and 
/o(980)7r+7r- 
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Fig. 1. Invariant mass distributions for K~ K + tt~tt + with standard cuts (a) used to 
determine the BR and tight cuts (b) used to study the subresonant structure. The 
fit (solid curve) is explained in the text, the dashed line shows the background 
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Fig. 2. Branching ratio T(D° — > ivT _ i ; r + 7r _ 7r + )/r(D — ► K _ 7r _ 7r + 7r + ) versus several 
sets of cuts. We varied the confidence level of the secondary vertex from 1% to 50% 
(16 points), Iso2 from 1CT 6 to 1 (7 points), L / <jl from 6 to 20 (15 points), Ak 
from 1 to 5 (9 points), picon from —6 to —2 (9 points), and finally we consider the 
tight cut set used to study the subresonant structure (last point). The dashed lines 
show the quoted branching ratio ±lcr. 
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ttV mass (GeV/c 2 ) K + "tt" + mass (GeV/c 2 ) 



Fig. 3. a) The K + K~ir~n + mass plot with tight cuts. The hatched area corresponds 
to the events used in the amplitude analysis. The other plots show two-body projec- 
tions of these events: b) K + K~; c) tt~tt + ; d) K + tt~ plus K~tt + (two independent 
entries per event). We observe clear signals of 0(1020), p(770)° and if* (892). 
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Fig. 4. The two-body projections with the fit result (solid line) and background 
(shaded histogram) superimposed: a) K + K~; b) 7r~7r + ; c) K + ir~ plus K~n + (two 
entries per event); d) K + tt + plus K~n~ (two entries per event). 
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Fig. 5. Two-dimensional projections: K + K~ vs. tt~tt + in the top row, K + tt~ vs. 
K~n + in the lower row. Data is in the left column; mini-MC, generated according 
to the fit solution, is in the right column. 
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Fig. 6. Distribution of cosine of the helicity and acoplanarity angles (see text for an- 
gle definitions) The bump in the central part in the distribution of cosine of p helicity 
angle (top right plot) is caused by the chain D° -» K 1 (1270)K, #i(1270) pK. 
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